Magnetization dynamics in the inertial regime: nutation predicted at short time scales 



O 

as: 
< 

(N 



M.-C. Ciornei and J.-E. Wegrowe 
Ecole Polytechnique, LSI, CNRS and CEA/DSM/IRAMIS, Palaiseau F-91128, France. 

J. M. Rubf 

Departement de Fisica Fonamental, Universitat de Barcelona, Diagonal 647, Barcelona 08028, Spain. 

(Dated: August 13, 2010) 

The dynamical equation of the magnetization has been reconsidered with enlarging the phase 
space of the ferromagnetic degrees of freedom to the angular momentum. The generalized Landau- 
Lifshitz-Gilbert equation that includes inertial terms, and the corresponding Fokker-Planck equa- 
tion, are then derived in the framework of mesoscopic non-equilibrium thermodynamics theory. A 
typical relaxation time r is introduced describing the relaxation of the magnetization acceleration 
from the inertial regime towards the precession regime defined by a constant Larmor frequency. 
For time scales larger than r, the usual Gilbert equation is recovered. For time scales below t, 
nutation and related inertial effects are predicted. The inertial regime offers new opportunities for 
the implementation of ultrafast magnetization switching in magnetic devices. 
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In 1935 Landau and Lifshitz proposed an equation for 
the dynamics of the magnetization M (of constant modu- 
lus) , composed of a precession term M x H and a longitu- 
dinal relaxation term M x (M x H) that drives the magne- 
tization towards equilibrium along the magnetic field H 
[1]. Two decades later T. L. Gilbert derived the equation 
that bears his name in which the relaxation towards equi- 
librium is described by a damping term -q [2] through the 
dynamic equation dWL/dt = 7M x (H — TjdM./dt), with 7 
as the gyromagnetic ratio. The two equations (Landau- 
Lifshitz and Gilbert) are mathematically equivalent. 

The range of validity of the Landau-Lifshitz-Gilbert 
(LLG) equation was established one decade later by W. 
F. Brown, with a description of a magnetic moment cou- 
pled to a heat bath ("thermal fluctuations of a single- 
domain particle", 1963 [3]). The magnetic moment is 
treated as a Brownian particle described by the slow de- 
grees of freedom (10~ 9 s), the angles {9, <p}. The remain- 
ing degrees of freedom of the system relax in a much 
shorter time scale (< 10~ 12 s). The time scale separa- 
tion between the rapidly relaxing environmental degrees 
of freedom and the slow magnetic degrees of freedom al- 
lows the coupling between the magnetization and the en- 
vironment to be reduced to a single phenomenological 
damping parameter 77, whatever the complexity of the 
microscopic relaxation involved [4, 5]. 

However, important experimental advances towards 
very short time-resolved response of the magnetization 
(sub-picoseconds resolution, i.e. below the limit proposed 
by Brown) have been reported in the last decade [6]. In 
parallel, industrial needs for very fast memory storage 
technologies are approaching the limits imposed by the 
precessional switching [7]. In these experiments and in 
the corresponding applications, time scale separation be- 
tween the conserved degrees of freedom {8, </>} and the 
other degrees of freedom, assumed by Brown [3] , finds its 
limit. 

The purpose of this Letter is to investigate the dynam- 



ics of the magnetization beyond this limit by extending 
the phase space to additional degrees of freedom expected 
to be also out-of-equilibrium at short time scales [5, 9]. 
According to the well-known gyromagnetic relation [10], 
the next relevant degree of freedom of the ferromagnetic 
system (beyond the coordinates of position; i.e. the an- 
gles {6, <fi}) is the angular momentum L. As will be 
shown below, the consequence of considering also the 
conservation of the angular momentum is that inertial 
terms, i.e. acceleration terms proportional to ePM/dt 2 , 
appear in the equation of motion. The existence of in- 
ertial terms in the dynamics of the magnetization opens 
the way to deterministic ultrafast magnetization switch- 
ing strategies, beyond the limitations of the precessional 
regime [11]. We assume however that the microscopic re- 
laxation channels are inactive at the time considered here 
(e.g. by choosing the adequate materials and excitations 
of the magnetization). Otherwise, a non-deterministic 
regime would take place [7, 12, 13]. 

We derive below the generalized Gilbert equation and 
the corresponding Fokker-Planck equation that includes 
the inertial effects for a uniform magnetic moment. The 
derivation is performed in the framework of mesoscopic 
nonequilibrium thermodynamics (MNET) [14-16], and is 
based on the expression of the conservation laws, ther- 
modynamic laws, and symmetry properties. 

It is convenient to model the dynamics of a magnetic 
moment M = M s e (submitted to an applied magnetic 
field H = — jr^^r an d coupled to a heat bath) with a 
statistical ensemble composed by non-interacting identi- 
cal uniform magnetic moments found in the same given 
conditions (ergodic property). Here, e, M s and V F are 
respectively the radial unit vector of angles {8, <f>}, the 
magnetization at saturation, and the ferromagnetic po- 
tential energy. The ensemble of magnetic moments of 
constant modulus M s defines a sphere surface £ and the 
number of magnetic moments oriented within (e, e + de) 
defines the density n(e) of magnetic moments over S. 
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We have shown in previous works that associating two 
degrees of freedom {9, (f>} to a magnetic moment is suf- 
ficient to derive both the Gilbert equation and the cor- 
responding rotational Fokkcr-Planck equation from non- 
equilibrium thermodynamics principles only [13, 17]. 

Extending the configuration space to the magnetic an- 
gular momentum L, the space £ is extended from a 
two dimensions space, to a priori five dimensions space 
{8, <fi, L} [18]. A distribution function /(e, L) of mag- 
netic moments with the magnetization orientation within 
(e, e + de) and the angular momentum within (L, L + c?L) 
should then be defined, where / is assumed to vanish for 
infinite values of L as: lim L ^±oo f( e , L) = 0. The an- 
gular momentum L associated to a magnetic moment is 
either changed by an applied torque N = M x H as 
(w) s = °i tner ^ ^e interaction with the heat bath. 
When considering the statistical ensemble, the interac- 
tion with the bath is modeled through a phase space flux 
Jl (defined below) which vanishes for large values of L: 
lim L ^ ±00 J L = 0. 

The kinetic energy expression that Gilbert associated 
to the magnetization [2] is written as: JC = LL : J _1 /2, 
where the magnetic inertial tensor /, is related to the 
magnetic moment (and not to the inertia of matter). It 
is assumed that / keeps the symmetry of the magnetic 
moment, i.e. is axial symmetric of symmetry axis e: I = 
Ii (l) — eej +I a ee, with U the dyadic unit (where I\ = I 2 
and I3 are the diagonal coefficients of the inertial tensor). 

In the space-fixed reference frame denoted by the sub- 
script s, the conservation law for the number of axial 
symmetric moments /(e, L) writes [19]: 
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where the derivatives with respect to the angles are made 
while holding the Cartesian components of L s constant: 



f d(fe) 
\ de 



sm( 



d(/sinflfl) 
06 



m, - 



The density n(e) of magnetic moments in the space X 
is recovered by integrating over the angular momentum 
degree of freedom n(e)= J /(e, L) d 3 L. The conservation 
law for the magnetic moments in the £ space is hence 
deduced from (1): 



dn_fdf_ 
~dt ~ J ~dt 



d 3 L s = -—- I fed 3 L s 



d(ne) 
de 



(3) 



magnetic angular momentum (L s ) is also derived [19, 20]: 



dt 



dt 

= J feL s d 3 L s +nN s (e) + J 



Jl d 3 L s 



d(L s ) d 



dt 



^- (exP s ) lnN s (e) + J ^d 3 h s (4) 



where the magnetic pressure tensor is defined as P s = 

T 1 J(L - (L a »(L - (L s ))/d 3 L s . 

The conservation equation (4) states that the rate of 
the average angular momentum (L s ) is due to three con- 
tributions: an applied torque N s , an average interaction 
with the bath /Jlc? 3 L s (i.e. damping), and a torque due 
to pressure (i.e. rotational diffusion). 

The expression for Jl is deduced from the entropy pro- 
duction expression cr(e) [14, 15, 20]. Defining the fer- 
romagnetic chemical potential ^(e, L), the power Ta(e) 
dissipated by the magnetic system is the product of the 
generalized flux by the generalized force: 



(5) 



where the chemical potential takes the canonical form 
[14, 16]: 

H(e, L) = kT In [f(L, e)] + /C(e, L) + V F (e) (6) 

The application of the second law of thermodynam- 
ics, together with the local equilibrium hypothesis in the 
(e, L) space, lead us to the introduce the Onsager matrix 
L such that: Jl = — C ■ As the Onsager coeffi- 

cients are a reflection of the system's symmetry [15], the 



relaxation tensor defined as r 



=-i 



jCI is also axial 



symmetric: r = t 1 _1 ([7 — ee) + r 3 ~ 1 ee (where n = r 2 
and T3 are the diagonal coefficients), and is related to 
damping. Moreover, as e is an axis of symmetry for the 
ferromagnetic potential V F (9,(j)), the relaxation tensor 

t is not expected to have any components in the e 
direction [19], leading to t^ 1 = 0. 

The dynamic equation (4) can be rewritten: 



d(L B ) 
dt 



As the inertial tensor I and the relaxation tensor r 
are time independent in the rotating frame (or magne- 
tization frame), a simpler expression of Eq. (7) can be 
obtained in this frame. After introducing the average an- 
gular velocity 17 such that (L) = I ■ CI, Eq. (7) rewrites 
as: 



Beyond, the conservation law for the mean value of the 



dtt r 
~~dT 



=-i 



n de 



P 



(8) 



3 



The rotating frame is denoted by the subscript r 
=-i =-i 



rot 



= T r 



Of- 1 ) 



1^3 e x U, or 



=-i 



(n a* 




and 



(9) 



where a* — a {h/h — 1) 1 with a — (Q3 n) 
The three components of Eq. (8) read: 
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. ^3 = — ' 3 

(10) 

Since the quantity L3 — is a constant of motion, 
the well-known gyromagnetic ratio 7 can be defined as 
the ratio of the magnetization at saturation M s by the 
axial angular momentum 7 = [2, 10]. 

Also, the averaged dynamic Eq. (10) introduces a char- 
acteristic time scale n, which separates the behavior of 
the magnetic system of particles in two regimes: the dif- 
fusion regime or the long time scale limit t » t\ , and the 
inertial regime or the short time scale limit i < ti. Since 
the modulus of the magnetization M is conserved, the 
relation ^ = fl x M holds. Cross-multiplying by M 
and using the above definition of 7 leads to the identity 

« = Mr X ' 



M 
" t " hi' 



M dM 

Mj A dt 

In a diffusive regime, i.e. for t ^> n , the inertial terms 
^ and ^p 2 - are negligible with respect to ^ and 
Eq. (8) then rewrites as the Gilbert equation with an 
inertial correction performed on the previously defined 
gyromagnetic coefficient 7* = — — 



1-/1//3 



dM 

— — = 7* M x 

dt ' 



H 



eff - V 



dM\ 
~dTj 



(11) 



The Gilbert damping coefficient r\ is now defined as: 
r] = ^ (so that a* = ^y*r]M s is the corresponding 

dimensionless coefficient), and H e // is an effective field 
that includes the diffusion term. 

At the diffusive limit, the magnetic moments follow a 
distribution function /(e, L) close to a Maxwellian cen- 
tered on the average angular momentum (L) [9]. This 
leads to a diagonal form for the pressure tensor: P = 
nkT/V and H e// = H - f ^ [13, 17]. 

Eq. (11) contains the density n(e) so that the equation 
is not closed. However, inserting Eq. (11) into the con- 
servation law (3) leads to the rotational Fokker-Planck 

equation of n(e), derived by Brown [3]: = a ^™^ n ^ 

For short enough time scales t w n , the inertial terms 
cannot be neglected and the Gilbert approximation is 
no longer valid. The dynamic equation (10) takes the 
following generalized form: 



dM 
~dT 
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7 d(M x P) 



n dM 

(12) 

The corresponding generalized rotational Fokker- 
Planck equation for the statistical distribution / is ob- 
tained with replacing J L by the Onsager relation derived 
earlier into the conservation law (1) and rewriting the 
law in the rotating frame [20]: 



dt 



= -i 

d(fe x I r L r ) 



+ 



d 

8L r 



fr ro \ ■ L r - /N r + kTr r ^r — 

(13) 



At short time scales t w n and due to the inertial 
effect, the usual precessional behavior is enriched by a 
nutation effect. The simplest way to understand nuta- 
tion is to imagine that the effective field is switched off 
suddenly with zero damping: the precession stops sud- 
denly because the Larmore frequency u>l = 7*H drops 
to zero at the same time. However, in the absence of 
inertial terms, the magnetic moment also stops at this 
position within an arbitrarily short time scale. But if the 
kinetic energy is not zero (and this is the case for magnc- 
tomechanical measurements of the magnetization [10]), 
the movement cannot be stopped suddenly: the preces- 
sion (around the magnetic field) stops but the magnetic 
moment starts to rotate around the angular momentum 
vector in order to conserve the energy: the precession is 
transformed into nutation. 

Fig. 1 shows the numerical resolution of Eq. (12) (ne- 
glecting thermal fluctuations) with a field along z axis 
and for a parameters n fixed to 2ps with \a\ = 0.05. 
The trajectories are plotted on the sphere E. The usual 
trajectory deduced from the LLG equation (n <§; ps) 
is also plotted for comparison. The motion of the mag- 
netic moment displays the familiar curve due to Larmor 
precession, with superimposed loops generated by the nu- 
tation effect. Fig. 1(b) shows a trajectory starting with- 
out initial velocity under an effective field of IMA/m, 
changed suddenly to 3MA/m and once again down to 
zero. Four curves are represented, two for the Eq. (12) 
with |a| = 0.05 (red continuous line) and a — (red 
dashed line), and two for the usual LLG equation with 
and without damping (blue). At the end of the motion 
(left), the field is set to zero and the precession is de- 
stroyed, with the nutation effect shaping a circle (with- 
out damping) or a spiral (with damping). Note that the 
profile of the nutation loops depends on the initial condi- 
tions (the cusp presented in Fig.l instead of loops is due 
to zero initial velocity). Fig. 1 (c) shows the time deriva- 
tive of the angle <j> as a function of time for the trajectory 
displayed in Fig. 1(b). The horizontal lines represent the 



4 




FIG. 1: Numerical resolution of Eq. (12) with n = 2ps. (a) 
Trajectories at two different fields with \a\ = 0.05 (red) and 
curves deduced from the LLG equation (blue), (b) Trajec- 
tory of the magnetization with changing suddenly the effective 
fields from H=l MA/m to H=3 MA/m and H=0, with damp- 
ing (continuous line) and without (dotted line), (c) Time 
derivative of the azimuth angle <f> plotted as a function of 
time for the trajectory of Fig. 1(b). 



constant Larmor frequencies, and the oscillations are due 



to nutation (for \a\ = 0.05 and a = 0). 

In the case of two stable states separated by a potential 
barrier (e.g. for magnetic memory units), efficient strate- 
gies based on the inertial mechanism can be perormcd for 
ultrafast magnetization reversal. Such a strategy has al- 
ready been implemented in the case of antiferromagnets 
[11], in which inertial terms are present due to the energy 
stored by deformation in the magnetic domain walls. The 
inertia has been used to overcome an energy barrier after 
having push the magnetization with a very short optical 
impulsion. The novelty of our results is that any kind 
of ferromagnets could in principle be used for ultra-short 
inertial magnetization switching. 

In conclusion, we have shown that extending the phase 
space of the magnetization to the degrees of freedom of 
the magnetic angular momentum leads to considere a 
generalized Landau-Lifshitz-Gilbert equation that con- 
tains inertial terms. This extension is justified by the 
well-known gyromagnetic relation that relates the mag- 
netization to the angular momentum. It is predicted that 
inertial effects should be observed at short enough time 
scales (typicaly below the picosecond), e.g. by measur- 
ing nutation loops superimposed to the usual precession 
motion of a magnetic moment. The inertial regime at 
short time scales would also offer possibilities for new ex- 
periments and devices based on ultrafast magnetization 
switching. 
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